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Abstract 

We examine the particle production via preheating at the end of inflation in super- 
symmetric theories. The inflaton and matter scalars are now necessarily complex 
fields, and their relevant interactions are restricted by holomorphy. In general 
this leads to major changes both in the inflaton dynamics and in the efficiency 
of the preheating process. In addition, supersymmetric models generically con¬ 
tain multiple isolated vacua, raising the possibility of non-thermal production of 
dangerous topological defects. Because of these effects, the success of leptogenesis 
or WIMPZILLA production via preheating depends much more sensitively on the 
detailed parameters in the inflaton sector than previously thought. 



1 Introduction. 


The inflationary paradigm [Q] has been remarkably successful in explaining the observed 
large-scale features of the universe. Apart from its prediction of a flat universe with 
VLtot = 1, inflation naturally produces primordial density fluctuations with a Harrison- 
Zeldovich power spectrum. These fluctuations in turn give rise to the angular pattern 
of the cosmic microwave background radiation (CMBR) consistent with recent observa¬ 
tions p. 

On the theoretical front, one of the most important recent developments in infla¬ 
tionary cosmology was the realization that non-perturbative particle production 

may play an important role in the inflaton decay process. The process of the inflaton 
decay through non-perturbative, parametrically enhanced particle production has been 
called preheating, to distinguish it from the usual reheating scenario where the inflaton 
decay is treated perturbatively. An interesting feature of the preheating scenario is that 
the production of particles with masses greater than the inflaton mass is possible. For 
example, it was claimed that in the simplest chaotic inflation model with the inflaton 
mass m ~ 10^^ GeV, the scalars with masses of up to 10^“^ — 10^^ GeVQ and fermions 
with masses of up to 10^^ — 10^® GeV can be copiously produced. In particular, the 
non-thermally produced particles could include the baryon number violating gauge and 
Higgs bosons of grand unified theories (GUTs) or right-handed neutrinos, resurrecting 
the possibility of GUT-scale baryogenesis or leptogenesis 0,^]. The non-thermally pro¬ 
duced superheavy particles could also play the role of dark matter p, and their decays 
could explain the observed super-GZK cosmic ray events 
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The magnitude of the observed GMBR anisotropies suggests that the energy scale at 
which inflation takes place is substantially lower than the Planck scale. Any microscopic 
model of inflation should explain this hierarchy of scales and its radiative stability. From 
theoretical point of view, the most well-established and appealing way to do this is to as¬ 
sume that the physics of inflation is supersymmetric. This assumption is independently 
motivated by the need to resolve the usual gauge hierarchy problem of the Standard 
Model, and by the consistency considerations in string theory. Supersymmetry puts 
non-trivial constraints on the spectrum of the theory: bosons and fermions necessarily 
come in pairs, and matching the bosonic and fermionic degrees of freedom is only pos¬ 
sible if all the scalar fields are complex. Moreover, the interaction terms possible in a 
supersymmetric theory are restricted by holomorphy. However, most of the studies of 
the preheating process to date concentrate on the case of a real scalar inflaton field, do 
not take into account the holomorphy constraints on its couplings, and consider pro¬ 
duction of fermions and bosons independently. In this paper, we would like to examine 
the inflaton dynamics and the process of preheating after the end of inflation, taking 
into account the constraints of supersymmetry. We will show that many interesting and 
important new features arise in supersymmetric models, which can substantially modify 


^Heavier scalar particles could be produced via a slightly different mechanism, “instant” preheat¬ 


ing 
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the conventional preheating scenarios. 

Let us mention that the parametric resonance for the case of a complex inflaton 
coupled to a scalar held was considered in Ref. |M. However, the structure of the scalar 


potential assumed in that study is quite different from the one studied here, which we 
believe is generic in supersymmetric theories. Also, we discuss preheating of fermions 
which was not considered in 
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The paper is organized as follows. We begin by reviewing the well-known qualitative 
results for preheating of scalar particles and fermions in the case of a real inhaton 
(which we refer to as “non-supersymmetric”) in Section We then proceed to discuss 
the general structure of the inhaton interactions in supersymmetric models, and explain 
the assumptions of our analysis (see Section |^.) In Sections § and ||, we discuss the 
preheating of scalar and fermion particles, respectively, in the supersymmetric case. We 
highlight the similarities and diherences between this case and the more familiar, non- 
supersymmetric case. We conclude in Section |^. Certain more technical points of our 
analysis are relegated to the Appendices. 


2 Review of Non-Supersymmetric Results. 

Before discussing the supersymmetric case, let us recall the standard theory of particle 
production via preheating in non-supersymmetric models. To be concrete, we consider 
the case of chaotic inhation with the inhaton potential R(0) = where the 

inhaton mass m = 10^^ GeV and the value of the inhaton immediately after 

the end of slow-roll inhation is 0o ~ Mp/3. We will review the non-thermal production 
of scalars and fermions in turn. 

2.1 Production of Scalar Particles. 

The production of scalar particles during preheating has been studied in detail in 0. 
Let us briehy summarize the results. Assume that the inhaton is coupled to a real scalar 
held X via 

= 2 (^0 + + 1 ^( 0 )- ( 1 ) 

In a hat FRW universe with a scale factor a{t), the equation for the mode of the held y 
with a comoving momentum k reads, 

Xk + 3Hxk + ^^kXk = 0, ( 2 ) 

where 

^ ( i -\2 + ^0 + ( 3 ) 

a(tY 

is the time-dependent ehective frequency of the mode and H = d/a. It is useful to 
hrst consider particle creation in Minkowsky space, with a = 1 and 0 = (po cos mt. The 
qualitative features of the process depend crucially on the parameter q = g'^(f)l/4:m‘^. 
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Typically g S> 1 during the preheating epoch. In this case, the y particles are created 
via broad parametric resonance. In this regime, the necessary condition for particle 
production is non-adiabaticity in the change of the effective frequency. 


du 

dt 




(4) 


This condition can be satished if 


+ e< 


mgcjyo 

2 


(5) 


That is, particles with bare masses of up to 


Mmax ~ \Jmg(j)o/2 = ( 6 ) 

are produced. For example, for g ~ 1, we hnd M^ax ~ 1-5 x 10^® GeV, more than two 
orders of magnitude above the inflaton mass. The maximum momentum for a particle of 
bare mass Mq is of order k^ax ~ Y^m(70o/2 — Mq . Note that for M close to Mmax, only 
particles with very low momenta can be produced, so that the production is suppressed; 
however, once M ^ Mmax/‘^, the typical momenta of the produced particles become of 
order Y^mg0o/2, and the number density of the produced particles becomes of the same 
order as in the case Mq = 0. 

The condition (^) is satisfied and the y particles are produced only for a short period 
of time during each oscillation, when the inflaton is close to the minimum of its potential: 


101 ^ 




< 00 - 


(7) 


For the rest of the oscillation period, the number of y particles remains constant. Because 
of this feature, the process of particle creation can be treated analytically |p. 

The expansion of the universe brings in two interesting new features. First, the 
number of particles with a given wavelength may both increase and decrease during a 
particular passage of the inflaton through its minimum, depending on the phase of Xk 
at that time. Nevertheless, the number of y particles is growing exponentially with 
time when averaged over many oscillation periods (this phenomenon has been called 
“stochastic” resonance 0) so preheating does occur. Secondly, the amplitude of the 
inflaton oscillations is damped with time due to Hubble friction, so that effectively the q 
parameter decreases with time. Particle production ceases when the condition (^) is no 
longer satished for any k, that is, when q = qmin ~ (dTo/m)^. (If Mq ^ m, the preheating 
stops when q ~ 1.) As we have mentioned above, for q ^ (2Mo/m)'^ = the rate 

of particle production is of the same order as in the case Mq = 0; for ~ g ~ ^Qqmin 
there is some production, but the rate is suppressed. 

In the above discussion, we have neglected the backreaction of the produced y par¬ 
ticles on the inhaton motion. If Mq = 0, this approximation breaks down when 


g ^ g2 


2giiMp j 

3m g^Mp 


( 8 ) 
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where /i is a parameter of order 0.1 which characterizes the efficiency of particle produc¬ 
tion. For most reasonable parameter values, the backreaction becomes important after 
0(10) inflaton oscillations [Q. The same estimate is true for massive particles, provided 
that their production is unsuppressed for all q > q 2 - This occurs if q 2 ^ lOgmm, or, 
assuming ~ 1, 

^0 < ~ 4m. (9) 

If, on the other hand, q 2 < qmin, the backreaction can always be neglected: this is the 
case for Mq ^ 8m = 8 x 10^^ GeV. Thus, there is a reasonably wide range of values of 
Mo (roughly between 10^"^ and 10^® GeV) where some particles are produced, but not 
enough for the backreaction to become an important issue. 

If the condition (D is satished, preheating proceeds in two stages. During the hrst 
stage, q > q 2 and the backreaction can be neglected. During the second stage, the 
backreaction has to be included. The “effective” g-parameter at this stage is given 
by geff = 0o/4(-^^)- This parameter is decreasing as more X particles are created: 
approximate energy conservation implies q^s oc ((X^))“^. Once q^s ~ qmin, particle 
production stops. 


2.2 Fermion Production. 


Several authors 
the Lagrangian 


12,13 


have studied the fermion production during preheating using 


Cf = (Mo g0)'0'0- 


( 10 ) 


They found that fermions as heavy as 10^^ — 10^® GeV could be produced. This is not 
surprising: with the Lagrangian as long as the inflaton oscillation amplitude 0o is 
greater than Mo/g, there is an instant of time during each oscillation when the effective 
mass of the fermion vanishes, and so the condition (^) is automatically violated at least 
for very large wavelengths. Thus, the upper bound on Mq is, for g ~ 1, of order of the 
initial value of the inflaton, 0.3Mp. 


3 Inflaton Interactions in Supersymmetric Models 

Now, let us discuss how the above discussion of preheating generalizes to the supersym¬ 
metric case. Interactions of the inflaton with other helds are of two types: those coming 
from the Kahler potential and the superpotential. While the Kahler potential couplings 
such as / d^eSX^X/Mp can play a role in perturbative reheating, they do not lead to 
resonant, non-perturbative particle production that we are interested in. Therefore, we 
will concentrate on superpotential interactions. The simplest supersymmetric general¬ 
ization of (0) is a model with two gauge-singlet chiral superhelds, the “inflaton” S and 
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( 11 ) 


the “matter field” X, whose interactions are described by the superpotentialQ 

W{S, X) = + 1(M + gS)X\ 

This is the most general renormalizable superpotential for two singlet fields with a 
discrete symmetry X —> —X. Without loss of generality, we will choose the parameters 
M, m and g to be real and positive. The scalar potential is given by 

V{S, X) = \mS + eS^ + ^gX^\^ + \M + gS\^\X\^, (12) 

where S and X are complex scalar fields. Note that this potential generically has four 
degenerate minima: {S = X = 0}, {S = —m/e,X = 0} and {S = —M/g,X = 
±^2{mM/g"^ — eM'^/g^)}. For e = 0, the number of minima is reduced to three: 

S' = X = 0 and {S = —M/g,X = ±\/2mM/g.} The presence of multiple, discon¬ 
nected degenerate vacua, which will have important consequences for the dynamics of 
preheating, is completely generic in supersymmetric theories with multiple singlets. In¬ 
deed, the scalar potential in such theories has a form 


r(*) = y 

i=l 


d(j)i 


(13) 


The vacua are found by solving a system of algebraic equations diWi^i) = 0,i = 
1.. .n. For a renormalizable superpotential, each of these equations is at most quadratic. 
Generically, the system has 2n different complex solutions, each of which corresponds to 
a vacuum with zero energy. One should keep in mind, however, that the degeneracy of 
these vacua is lifted by supersymmetry breaking unless the theory possesses additional 
discrete symmetries relating different vacua. (The model defined by (|^) happens to 
have such a symmetry, the X —>■ —X reflection.) 

The expression (0) can be thought of as the leading terms in the Taylor expansion 
of the full infiaton potential around S' = 0. This potential is sufficient to describe the 
(p)reheating phase considered here, and is independent of the details of the slow roll 
phase of inflation.^ The only assumptions we make concern the values of the S and X 
fields at the end of the slow roll phase, which set the initial conditions for the process 
we are considering. Specifically, we assume that the initial value of the infiaton field S'o 
is somewhat lower, but not too far from, the Planck scale Mp, while the initial value of 
X is zero. The latter assumption is reasonable since the mass of X is large (Planckian) 
during inflation. We treat the initial phase of the complex field S as arbitrary, which is 

^Supersymmetric models of chaotic inflation were constructed in Ref. [ p^ . 

^We expect the qualitative features of our analysis to hold in both chaotic inflation and new infla¬ 
tion models. Parametric resonance in supersymmetric hybrid inflation models was considered in [ p^ ; 
however, it was subsequently shown that in these models the infiaton typically decays via tachyonic 
preheating |j^, without developing a resonance. 
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a good approximation as long as the potential in the phase direction is sufficiently flat 
during slow-roll inflation. We will study the dependence of the results on this phase. 

The remaining freedom concerns the choice of parameters in the superpotential (pd^. 
Throughout the analysis, we will assume that m <C l^ol; for numerical estimates, we 
will use m ~ 10^^ GeV and l^ol ~ Mp/3 to facilitate the comparison with the non- 
supersymmetric analyses reviewed in Sec. (These particular choices are motivated by 
the models of chaotic inflation [|l|.) In this paper, we will primarily be interested in the 
production of particles that are too heavy to be produced through thermal reheating or 
perturbative inflaton decay; therefore, we will assume M > m. The size of parameter 
e in chaotic inflation models is constrained to be at most about 10“® lITTll. Even such 


small coupling, however, can have profound effects on the inflaton evolution during the 
(p)reheating epoch. We will therefore consider two qualitatively different, physically 
interesting cases: e = 0 and e ~ 10“®. In the first case, the complex phase of the 
inflaton is conserved during the evolutionj^ and the motion of the inflaton is identical 
to the non-supersymmetric case. Still, there are some important differences as far as 
production of particles (especially fermions) is concerned. In the second case, the inflaton 
generally acquires angular momentum and moves along an elliptical spiral trajectory in 
the complex plane. As a result, the preheating process is quite different from the non- 
supersymmetric case. 


4 Scalar Particle Production in Supersymmetric Mod¬ 
els. 

In this section, we will consider the resonant production of scalar X particles during the 
reheating era in the model defined by the superpotential (pT|). 

4.1 e = 0. 

We start by examining the situation where the coefficient of the S'^ term in the super¬ 
potential (^) vanishes. As long as the backreaction of the X particles can be neglected, 
the motion of the inflaton is described by the potential V{S\,S 2 ) = m?\S\‘^. Since at 
the end of inflation the inflaton is slowly rolling, i.e. S' ~ 0, its subsequent evolution is 
given by 

S{t) = S'oe*^ cosmt, (14) 

where 6^ is a constant phase, determined by the initial conditions which we will assume 
to be completely random. In an expanding universe, the amplitude of the oscillations 
S'o becomes time-dependent; assuming S'o ~ Mp/3 at t = 0, after the first few oscil¬ 
lations So ~ Mp/3mt ^ Mp/20N, where N is the number of oscillations that have 

"‘Strictly speaking, this is only true as long as the backreaction of the produced X particles can be 
neglected - see Sec. 
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been completed. The motion of the inflaton described by ([T^) is identical to the non- 
super symmetric case, see Sec. 

The mass terms for the X held can be read from (1^): 


^(M -|- gSiY + 9^82 + gmSi^ + (^{M -|- gSi)^ + g'^S^ — gmS^ X2 -h 2mgS2XiX2, 

(15) 

where we have decomposed X = Xi + 1 X 2 and S = Si + 182 - Diagonalizing this mass 
matrix gives 

((M + gSi)^ + g^Sl + Xl + ((M + gSif + g^Sl - Xl (16) 


where Xi ^2 are the mass eigenstates. (In what follows, we will drop the tilde to avoid 
cluttering.) The mode equations for the helds Xi 2 have the same form as (H), with the 
effective frequencies 

981 ^ + ± gm\S\, (17) 

where the upper (lower) sign applies to the mode Xi {X 2 ). This equation has some 
important differences from its non-supersymmetric counterpart, Eq. (|^). In particular, 
the right-hand side of ([l7|) is not positive-dehnite. Since we are principally interested 
in heavy particle production, let us assume that m -C M. Then, the right-hand side 
of (prl) for the held X 2 is always positive, whereas for the held Xi it can become negative 
for sufficiently small k. The region in the Si — 82 plane in which this is possible (the 
“instability region”) is given by 


{Si - S,f + Si < 


mM 

9 ’ 


(18) 


where S'* = —M/g^ and terms of order m/M have been dropped. The instability region 
and a typical inhaton trajectory in the e = 0 case are shown in Fig. 1. 

Just like in the non-supersymmetric case, resonant production of X bosons occurs 
whenever the condition (|^) is violated. For the mode X 2 , this can happen if 


TTl 1 

M^(sin^ ^ M ~ 29 ^ 80 . 


(19) 


For the mode Xi, the analysis depends on the phase 6 . If 9 ^ iJm/M, the condition is 
similar to ( 0 ): 

777 1 

M^(sin^ 6^ — — cos^) ^-^fmS'o. (20) 

If, on the other hand, 9 < the inflaton will pass through the instability re¬ 

gion (|I8|), provided that the oscillation amplitude is large enough, S'o > IS"*!. Dur¬ 
ing this passage, the adiabaticity condition is bound to be violated, and production of 
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Figure 1: Motion of the inflaton in the Si — S 2 plane in the case e = 0. The red circle 
indicates the instability region (^). 


Xi particles will occur. Thus, scalar particles with bare masses as high as qSq (i.e. 
up to about 10^® GeV, assuming g ^ 1) can in principle be produced in our model. 
However, production of particles substantially heavier than the limit indicated by the 
non-supersymmetric analysis of Sec. ||, Eq. (|]), is only possible if 0 < Since 

m -C M, the range of 6 for which this is true is quite small; in other words, production 
of superheavy scalar particles requires fine-tuned initial conditions. 

The resonant production of X particles can be described analytically using the tech¬ 
nique developed in Ref. |^. The relevant calculations are presented in Appendix |^. 
For generic initial conditions on the inflaton, the production rates are identical (up 
to 0{m/M) corrections) to those found in the non-supersymmetric case with the po¬ 
tential (|l|) and Mq = M sin 9. For special initial conditions, when the inflaton passes 
through the instability region, the production of Xi particles is enhanced due to their 
locally imaginary mass. As long as the inflaton oscillation amplitude is large, Sq 3> |S'*|, 
the Xi production rate is close to the rate estimated in for non-supersymmetric par¬ 
ticles with zero bare mass (see Eq. with Mq = 0.) On the other hand, when Aq ~ S'*, 
the Xi production rate can be dramatically enhanced.^ These results are derived in 
Appendix 

To summarize our discussion so far, for generic initial conditions on the inflaton field, 

^Let us comment on the relation of our analysis in the case of “special initial conditions” to the 
tachyonic preheating scenario considered in [^. In both cases, the field coupled to the inflaton develops 
an instability. However, in our case, the instability occurs only for a very short time period during each 
inflaton oscillation. As a result, the fractional energy loss by the homogeneous component of the inflaton 
in each passage is small. (This is true as long as Sq ^ 5”* - a condition that is typically satisfied for 
a large number of inflaton oscillations at the beginning of the reheating era.) The backreaction of the 
produced X particles does not become important until at least several inflaton oscillations have taken 
place, and the language of stochastic resonance is applicable. 
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0 ^ Jm/M, both Xi and X 2 particles are resonantly produced provided that 


M < 



( 21 ) 


This condition is similar to the one obtained in for the non-supersymmetric case, 
Eq. (H). Indeed, if (U) is satished, preheating occurs regardless of the initial phase of the 
inflaton held. Preheating of particles with higher bare masses is possible, but requires 
special initial conditions. The absolute upper bound on M is achieved for the case when 
the inhaton is real {9 = 0) and is of order qSq ~ 10^^ GeV for ~ 1. Moreover, for 
the same special initial conditions, the production of Xi bosons can be signihcantly 
enhanced. 

The above discussion neglected the backreaction of the produced X particles on 
the motion of the inhaton. Let us now include this ehect. First, consider generic 
initial conditions, 6 S> \Jm/M. In this case, the rough estimate of the produced X 
number density is the same as in the non-supersymmetric case discussed in Sec. H, 
with the replacement Mq —> M sin 9. Thus, the backreaction can be neglected for 
M sin 9 ^ 8m 8 x 10^^ GeV. For smaller values of Msin9, the X particles are still 
being produced at the time when their backreaction becomes important. Including it 
leads to an ehective inhaton potential of the form 

K(S) = / (|Xp) ((Si - S. f + si) + m2|S|2. (22) 

The phase of S is not conserved by the hrst term in this potential, and will start changing 
(tending towards zero) when the backreaction becomes important. As a result. Si and 
S 2 oscillations will be out of phase. The resonant X particle production will quickly 
terminate, since it is possible only when Si and S 2 are both close to the minimum of 
the potential. For the modes Xi 2 (A:) with large occupation numbers, the subsequent 
evolution can be considered classically. However, it is still very complicated, with all the 
four coupled helds Xi 2 and 51^2 undergoing large oscillations in the potential 


V{S,X) « ((Si - S.)" + S|) + \g^X\\ (23) 

where we have assumed (|Xp) 3> wS. The oscillations are damped by Hubble fric¬ 
tion. When the oscillation amplitudes become sufficiently small, (|Xp) < rrS' or 
{{Si — < mM, the terms in the scalar potential which were neglected in Eq. (^31) 

become important. It is at this stage that the system chooses which of the three vacua 
will be the hnal point of the evolution. This is a very important question. If the vacua 
with non-zero X and S are preferred, the phenomenology will crucially depend on the 
exact nature of the X held. For example, if the X held carries R parity, these vacua 
are clearly undesirable. Moreover, in our model there are two vacua with non-zero X. 
The X helds whose evolution is considered here are position-dependent, since modes 
in a wide range of wavelengths obtain large occupation numbers during preheating. It 
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is therefore very likely that different vacua will be chosen in different spatial regions, 
leading to formation of domain walls. (This non-thermal defect production mechanism 
is similar to the one considered, for example, in Because of the Z 2 symmetry relat¬ 
ing the two vacua, they will remain degenerate even after supersymmetry is broken, and 
the domain walls of our model are stable.|^ Therefore, the model is only consistent with 
standard cosmology if the system ends up in the vacuum S' = X = 0 throughout the 
space. It would be very interesting to study the evolution in more detail and understand 
which vacuum is chosen for various values of the parameters and initial conditions. This 
will likely require a numerical study using the methods developed in 

For special initial conditions, d ^ the backreaction is expected to be very 

important. Indeed, the Xi particle production rate in the instability region (p!8D is at 
least as high as the rate in the non-supersymmetric, Mq = 0 case considered in 0. 
Based on the analysis of |^, we conclude that the backreaction of Xi particles will be¬ 
come important after just a few passages of the inflaton through the instability region. 
Moreover, in this case the Xi production will continue even after the backreaction be¬ 
comes important, and only stop when (|Xp) ~ {{Si — 5*)^). After that, the evolution 
of the modes Xi{k) with large occupation numbers can be described as classical motion 
of the helds Xi and Si in the potential (p3|), with X 2 ~ S '2 ~ 0, damped by Hubble 
friction. A quick estimate suggests that the evolution is more likely to end in the vacua 
with non-zero X, leading to domain wall production. 


Summarizing, we have found that for generic initial conditions, 6 S> ym/M, and 
sufficiently high masses, M sin 6 ^ 8 m, preheating shuts down before the backreaction 
of the produced particles on the inflaton motion becomes important, and the picture 
is quite similar to the non-supersymmetric case studied in [^. On the other hand, 
for MsinO ^ 8m the backreaction becomes an important effect. Once it is taken 
into account, the evolution of the system is rather more involved than in the non- 
supersymmetric case, and may lead to unacceptable consequences such as production 
of domain walls. Clearly, further investigation of this situation is necessary. Further¬ 
more, for special initial conditions, 9 ^ ^m/M, we found that the upper bound on 
the masses of the X bosons that can be reheated is substantially higher than in the 
non-supersymmetric case. In this case, backreaction is expected to be important, and it 
seems likely that domain walls are produced at the hnal stages of the process. Again, a 
more detailed numerical study is necessary to check this conclusion. 


4.2 e ^ 0. 

The cubic term S^ in the superpotential (pT]) has important consequences for the motion 
of the inflaton after the end of the slow-roll inflation and, therefore, for preheating. In 
the presence of this term, the motion of the inflaton is governed (neglecting backreaction) 

®In more general models, supersymmetry breaking can lift the vacuum degeneracy, leading to unsta¬ 
ble domain walls. In this case, the cosmological bounds could be avoided, but a more careful examination 
of the domain wall dynamics is necessary. 
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by the potential 


2\2 


1/(^1, ^ 2 ) = m\Si + ^ 2 ) + 2em(5^ + ^ 2 )*^! + + ^ 2 ) 


(24) 


The second term in this potential violates the symmetry S' —>■ e^^S. Due to this term, 


the motion of the inflaton for generic initial condition is no longer decribed by 
but is much more complicated. We will restrict our analysis to the case when the terms 
proportional to e in (1^) are of the same order as the hrst term at the end of the slow-roll 
inflation. (With our standard assumption m ~ 10^^ GeV this implies e ~ 10“®, close to 
its upper bound in chaotic inflation models [|^.) The evolution of the inflaton during the 


reheating era can be broken into two stages. During the hrst stage, the effect of the cubic 
and quartic terms in the potential is signihcant, and the motion is complicated. The 
potential (|^) has two vacua: S' = 0 and S = —m/e. Because of the Hubble friction, the 
energy of the system is decreasing and it will eventually settle down into one of these 
vacua. Let us assume that the system ends up in the vacuum with S' = 0. (This is 
always the case when m/e 3> |S'o|. For the case considered here, m/e ~ |S'o|, the choice 
of vacuum depends on the initial conditions; domains of attraction of the two vacua 
have more or less the same size.) In this case, the average magnitude of the inflaton held 
I S'! = \JSi + S 2 decreases with time. Eventually the quadratic term in (^) becomes 
dominant, and the second stage of the inhaton evolution begins. At this stage, the 
potential can be approximated by setting e = 0. However, the problem does not reduce 
to the case studied in the previous subsection, since the inhaton held generically has 
non-zero angular momentum. Its motion can be described as independent, out-of-phase 
oscillations of Si and S 2 helds: 


*S'i = S'i,o(^) sinmf, 

S 2 = S 2 p(t) sm{mt + S), (25) 

where and S' 2,0 are slowly decreasing with time as a result of the Hubble expansion. 
Thus, during the second stage, the inhaton is approaching the minimum of its potential 
along an elliptical spiral. As we show in Appendix P, with the assumptions made here, 
the ellipticity of the orbit is generally of order one. A typical trajectory of the inhaton 
is sketched in Fig. 2. 

The resonant production of X particles is only possible if the non-adiabaticity con¬ 
dition (^ is satished. For low-momentum modes, the frequency can be approximated 
by 

= {MgSifg'^S^T gmlSl- (26) 

During the hrst stage of the reheating era, the resonant X production is very unlikely. 
Indeed, approximate energy conservation implies that the velocity of the inhaton satishes 
\S\ ^ mlAol. Since m <C fi'|S'o|, the non-adiabaticity condition can be satished only if the 
frequency oj is much smaller than its typical value, which is order g\So\. If M g\So\, 
this can only occur when both Si and S '2 are close to zero, |S'| ^ M; if M 5 f|S'o|, this 
requires that the inhaton pass through (or close to) the instability region (|T^. The area 
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Figure 2: A typical trajectory of the inflaton in the case 
trajectory depends on the initial conditions and does not 
second part of the trajectory, corresponding to els'! m, 
the phase of S vanishes at the end of inflation.) 


e 7 ^ 0. The hrst part of the 
have a universal shape. The 
is an elliptical spiral (unless 
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Figure 3: Depending on the step of the spiral, the inflaton either passes through the 
instability region (left panel) or misses it (right panel). 


of these “production regions” in the Si — S 2 plane is quite small compared to the range 
of the inflaton motion during the hrst stage. Since the inflaton motion is highly non¬ 
periodic, it is very unlikely that it will pass through the “production regions” repeatedly, 
as required for the buildup of the number density in the resonant production scenario. 

Now, let us consider the second stage of the preheating era, when the motion of the 
inflaton is approximately described by (|^). For simplicity, let us set S' 1^0 = S 2 fl and 
6 = 7r/2, so that the orbit is a circular spiral. (The physical results of our analysis 
are insensitive to these assumptions.) In this case, the condition (H) can be studied 
analytically. The bottomline is that the X particle production only occurs when the 
inflaton passes through the “production region” in the Si — S 2 plane. This production 
region coincides (up to order-one factors) with the instability region dehned by Eq. (pIS]). 
This is not surprising: it is clear that the adiabaticity condition is always violated in the 
instability region. In the e = 0 case studied in the previous subsection, we found that 
X particles can be produced even without the inflaton passing through the instability 
region as long as their mass is sufficiently low, see Eq. (^). This production occurs 
twice per oscillation period when the inflaton passes the origin, provided that its velocity 
1^1 mSo is greater than (roughly) M/g. In the case of spiral motion considered here, 
however, the inflaton does not get close to the origin until it loses most of its energy. 
Thus, the only possible production mechanism in this case is by passing through the 
instability region. 

The inflaton will pass through the instability region regardless of the initial conditions 
only if the step of its spiral trajectory 5 at the time when IS"! ~ S'* is much smaller than 
the extent of the instability region along the Si axis, AS ~ yrnM / g. (This is illustrated 
in Fig. 3.) Estimating 5 ~ S^/Mp, we hnd that this requirement implies an upper bound 
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on the bare mass of the X particles that can be prodnced: 


M ;£ 


( 27 ) 


This bound is somewhat higher than the corresponding bound in the non-supersymmetric 
case, Eq. (|^), and the bound for generic initial conditions in the case e = 0, Eq. O' 
For our standard parameter values, ~ 1 and m rsj 10^^ GeV, Eq. (pTl) implies M < 10^^ 
GeV. 

Particle production by a spiralling inflaton passing through the instability region is 
considered quantitatively in Appendix^. The production is extremely efficient: in fact, 
the backreaction of the produced Xi particles on the inflaton may become important 
already after two or three production events, terminating the process. This explosive 
production is analogous to the so-called tachyonic preheating IM, which frequently 


occurs in hybrid inflation models. The subsequent evolution of the system can be studied 
classically, but is very complicated due to a large number of helds involved. It seems 
likely that this evolution will involve domain wall formation. 

In summary, we have found that if the cubic term is present in the inflaton superpo¬ 
tential, the dynamics of preheating is qualitatively different from the non-supersymmetric 
case. If the condition (|^) is satished, non-perturbative particle production occurs when 
the inflaton trajectory is spiral, and does not require hne-tuned initial conditions. Pro¬ 
duction in this case tends to be explosive, with the inflaton losing a large part of its 
energy already during the hrst two or three production events. For larger values of 
M, particle production may in principle occur during the initial, irregular phase of the 
inflaton motion; however, for generic initial conditions, this possibility is highly unlikely. 


5 Fermion Production in Supersymmetric Models. 


Preheating of fermionic particles in supersymmetric models has important differences 
from the non-supersymmetric case reviewed in Sec. |2.2|. The matter chiral superheld X 


contains two two-component (Weyl) spinors, which we will denote by 'ifji and 'ifjR. The 
interactions of these helds follow from the superpotential (pdj) : 


Cf = -(M + - (M + 


(28) 


Gombining the helds 'ipi and i/jr into a four-component (Dirac) spinor T = ■, 

Eq. (pSD can be rewritten as 

Cf =-{M + gSi)^^ + (29) 

where S' = S'! -|- zS 2 , T = and we have chosen the Weyl basis in which 7 ® = 

diag(—1, —1,1,1). Gomparing Eq. (pUD with its non-supersymmetric counterpart, Eq. O). 
indicates an important diherence: the appearance of an additional, pseudoscalar mass 
term in the supersymmetric case. This term plays a signihcant role in the analysis of 
preheating. 
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The evolution of the fermion wavefunctions is governed by the Dirac equation of the 
form =A(f)T(t), where 

±/c M + gS{t) \ 

M + gS*{t) Tk )' 

In this equation, k is the momentum of the mode and we are working in the helicity basis, 
ia^diip = ilk'll!. The frequencies of the fermionic modes are given by the eigenvalues of 
the matrix A, which are in general time-dependent: ui = i\Jk‘^ + \M -|- gS{t)\‘^. A 
necessary condition for the adiabatic evolution of the system is ci; -C o;^; when this 
condition is violated, the evolution is non-adiabatic and particle production is expected 
to occur. Below, we will use this condition to obtain a qualitative picture of fermion 
preheating in supersymmetric models. 




5.1 e = 0 


As expected in a supersymmetric theory, the analysis of fermionic preheating in this case 
is very similar to the scalar preheating analysis of Sec. B For m -C M < ^ mg So/2, 

fermions can be produced for any value of the inflaton phase, whereas for mg So/2 ^ 

M ^ gSo, the production occurs if this phase satishes 


sin^l ^ 


1 j gmSo 

MV 2 


(31) 


The production of very heavy fermions (up to M ~ gSo) is possible,!] but only if the 
initial phase of the inflaton is hne-tuned to be close to 0. (This hne-tuned case is identical 
to the non-supersymmetric model of Sec. |2]^.) For generic initial conditions, the upper 
bound on the mass of the fermions that can be effectively preheated is of order ^JgmSo/2, 
well below the non-supersymmetric estimate. 

The condition D <C o;^ is necessary for adiabaticity, but it is not sufficient. Particle 
production could occur even when D = 0 if the eigenvectors of the matrix A have strong 
time-dependence. For example, a rapid change of sign of the effective mass M -|- gS, 
even without a change in its magnitude, would lead to particle production. In the case 
e = 0, however, it is possible to show that the eigenvector rotation is adiabatic whenever 
the condition ci; -C o;^ is satished. Thus, fermion production does not occur outside of 
the parameter regions discussed in the previous paragraph. 

It is important to keep in mind that in supersymmetric models, fermion and boson 
production cannot be considred separately: the fermion and scalar masses and their cou¬ 
plings to the inflaton are all expressed in terms of just three superpotential parameters, 
M, m and g. As we emphasized in Section f4.1| , some regions of this parameter space 

^Naively, the expression (H seems to imply that infinitely heavy fermions can be produced for 
0 = 0. This is of course incorrect; the bound (0) is not applicable for 9 < ^m/M. The correct bound 
in that case is M < gSo. 
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may not lead to consistent cosmology, for example due to non-thermal production of 
topological defects. If this is the case, additional restrictions would be imposed on the 
masses of the fermions that can be preheated. In particular, the production of super¬ 
heavy fermions, M > y/gmS^, requires a passage of the inflaton through (or close to) 
the instability region, which results in enhanced scalar particle production and may be 
accompanied by domain wall formation. 

Even if the fermion is light enough to be preheated, we expect the number of produced 
particles for generic initial conditions to be smaller than in the non-supersymmetric case. 
Indeed, the effective frequency of the modes in the supersymmetric case is limited from 
below by M sin 9, leading to a suppression in the efficiency of particle production for non¬ 
zero 9. The magnitude of this suppression could be estimated by a numerical integration 
of Eq. (|30[) . Such an investigation is outside the scope of this paper. 


5.2 e y 0. 


Again, the analysis is completely analogous to the corresponding analysis of the scalar 
production presented in Section |4.2| . There are two possibilities for fermion production. 
First, it can occur when the inflaton passes close to the origin in the Si — S 2 plane, with 
a sufficiently large velocity. The second possibility is when the inflaton passes through 
the “instability region” around the point Si = —M/g,S 2 = 0. (The fermion effective 
mass vanishes at that point.) To obtain resonantly enhanced production, the inflaton 
should pass through either one of these special regions repeatedly. 

For non-zero e, the reheating epoch can be subdivided into two stages. During the 
hrst stage, the inflaton motion is highly non-periodic and irregular (see Fig. 2.) There¬ 
fore, repeated passage through either one of the two relatively small production regions 
is unlikely. During the second stage, the inflaton approaches its minimum along an ellip- 
tically spiral trajectory. Fermion production is possible if this trajectory passes through 
the instability region, given by (pISl) . As shown in Section O, this puts an upper bound 
on the bare mass of the particles that can be preheated, M < . However, if 

this condition is satished, the production of scalars in the instability region is explosive, 
generically leading to undesireable consequences such as domain wall formation. Thus, 
obtaining a cosmologically consistent scenario of fermionic preheating in a supersymmet¬ 
ric model with e 7 ^ 0 is difficult. Given this conclusion, we will not attempt to calculate 
the number density of fermions produced in this case. 


6 Conclusions. 

In this paper, we have studied the implications of supersymmetry for the inflaton dy¬ 
namics in the reheating era and the process of preheating. The emerging picture shows 
important differences from the previous, non-supersymmetric analyses. This is due to 
several robust, generic features of supersymmetric theories. First, all the scalars in su¬ 
persymmetric theories are necessarily complex, meaning that an extra degree of freedom. 
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the inflaton phase, has to be considered. Second, the conplings between the inflaton and 
the “matter” fields (the fields that are being preheated) are constrained by holomor- 
phy. Third, snpersymmetric models nsnally possess mnltiple degenerate vacna, which 
can become relevant (and lead to cosmological difficnlties) if the density of particles 
prodnced dnring preheating is large. Finally, the preheating of bosons and fermions in 
a snpersymmetric model cannot be considered in isolation: both kinds of particles are 
necessarily present, and their conplings are related. 

Several applications of preheating were snggested in the literatnre, and have been 
argned (nsing non-snpersymmetric models) to be snccessfnl. Do these optimistic con- 
clnsions snrvive in snpersymmetric theories? Onr analysis demonstrates that the answer 
to this qnestion is not straightforward. For example, consider a scenario of leptogene- 
sis where heavy (10^^ GeV) right-handed nentrinos are prodnced via preheating. This 
process is described by onr toy model, Eq. (^), if we identify X with the right-handed 
nentrino chiral snperfield. Generically, the parameter e is non-zero, and the analysis of 
Sections 4.2 and 5.2 applies. Eq. (|27|) shows that snfficiently heavy right-handed nen¬ 
trinos can indeed be prodnced. However, at the same time, their scalar snperpartners 
are prodnced explosively (see Section 4.2 and Appendix A.2.) As a resnlt, the system is 
likely to end np in the vacna with non-zero right-handed snentrino vacnnm expectation 
valne, either in some regions of space or everywhere. Both possibilities are phenomeno¬ 
logically disastrous. Therefore, the overall success of the leptogenesis scenario is not 
guaranteed, and may involve significantly more parameter fine-tuning and/or additional 
model building than previously thought. 

Our analysis is only the first step in the quantitative investigation of preheating after 
inflation in supersymmetric models. It needs to be extended in several directions. For 
example, in many cases, preheating is so efficient that the backreaction of matter should 
be taken into account at the end of this process. We have neglected this effect throughout 
the analysis, giving only qualitative estimates for the parameter ranges over which it 
becomes important. Studying the evolution of the inflaton-matter system after the end 
of preheating, including the backreaction, is crucial since it will allow one to answer the 
all-important questions about the vacuum selection and domain wall formation. Also, 
we have not attempted to describe quantitatively the regions of parameter space that 
allow for successful applications of preheating, e.g. to leptogenesis or WIMPZILLA 
production. Further work in this direction is necessary. 
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A Analytic Description of Scalar Preheating. 

In this appendix, we will present a semi-analytic treatment of the production of scalar 
particles X during inflaton oscillations. The qualitative aspects of this phenomenon have 
been discussed in Section |[ The calculation presented here is very similar to the one 
performed in P] for the case of real scalars. However, our results are more general; in 
particular we will consider the case when the inflaton field passes through the instability 
region (^), and show that this generally leads to enhanced particle production. 


A.l e = 0. 


Let us hrst concentrate on the case e = 0, discussed in Section ^?T|. Consider production 


of the quanta of the field Xi. Defining 0(f) = a ^/^(f)Xi(f), the mode equation reads 


with the effective frequency 

k2 


= 


a^(t) 


©fc + — 0 


g\Si - S^y + - gm\s\ + A, 


(32) 


(33) 


where S'* = —M/g, and terms of order m/M are neglected. The correction term A = 
— \{a/ay — I (a/a) is also typically very small and will be neglected throughout the 
analysis. 

The motion of the inflaton held for e = 0 is described by Eq. (|T^. Particle production 
occurs only when the non-adiabaticity condition (^ is satished. For generic initial 
conditions, this requires M <C Sq, see Eq. and production occurs twice during each 
inhaton oscillation cycle, when S is close to the origin. For special initial conditions 


production occurs when the inhaton passes through the instability region ([T8|) . In both 
cases, each of the production periods is short compared to the inhaton oscillation period. 
This leads to the following approximation. Let tj denote the moments of time when the 
S'-dependent terms in are minimized: 

McosS 


COS mtj = 


gSo 


Particle production occurs at tj — At ^ t ^ tj + At, where mAt 1. Outside of 
these time intervals, the mode equation (l3^) can be solved by making the adiabatic 
approximation. For fj_i + At ^ t ^ tj — At, we write 


Qk{t) = 


(y.1 




exp ( —i / ojdt ) -1- 




\p2uj 


exp ( / ojdt 


(35) 


where the coefficients a], and Pi are constant. In the next adiabatic time period, tj+At 
t ^ fj+i — At, the expression for 0fc(f) has the form 


< 

rsj 


Qk{t) = 


a 


i+i 


■\/2uj 


exp 


-i J ujdt 


+ 




:i+i 

k 


\p2uj 


exp ( +i / oodt 


(36) 
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With our normalization, the coefficients a and (3 are just the coefficients of the Bo- 
golyubov transformation of the creation and annihilation operators. The comoving den¬ 
sity of Xi particles at time t is given by 

1 


When the adiabaticity condition holds, the number density in (|T7D is constant, since Pk 
does not change. When the condition is violated, the number density can change: in 
general 7 ^ PP Thus, to estimate the change in particle density during a single pro¬ 
duction event, we need to hnd the Bogolyubov transformation matrix B corresponding 
to this event: 



(38) 


To achieve this, we perform a Taylor expansion of the effective frequency (|33|) for t k, tj\ 


00 


2 

k 






(39) 


where ~ sirP 9 — mM cos 9. For generic initial conditions, the production 

occurs when S' ~ 0, and therefore |*F(fj)| ~ mS'o; whereas for special initial conditions, 

we have |^(tj)| ~ Sq — S^. Dehning fc* = \Jg\S{tj)\, r = kpt — tj), n = k/a{tj)k^, 
and using the mode equation at t ~ tj becomes 


d^Qk 

dr'^ 


+ (H + 


0 


(40) 


where 

d ~ cos 9^ . 


(41) 


Note that the parameter XI can be negative for modes with low momenta if the 
passes through the instability region (|^. The exact analytic solution of (^OD 
by parabolic cylinder functions m-- 


inflaton 
is given 


0fc(r) = + z)t) + C^Dpi-a + *)r), (42) 

where p = —(1 -|- iX^P)/2 and Ci ^2 are arbitrary coefficients. The form (^21) is valid 
for both positive and negative A. The asymptotic forms of 0fc('r) for r —>■ —00 and 
r —>• -|-cx) should match the adiabatic solutions before and after the production event, 
Eqs. (|35|) and (|36D respectively. Performing this matching allows us to express the four 
Bogolyubov coefficients, and in terms of only two constants Ci and C 2 , and 

therefore to hnd the transformation matrix B;(, in (^). The result is 


( Vl + j^-^x/2+2iei 


(43) 
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where A = X{, 0{ = Jq uj{t)dt is the phase accumulated by the moment tj, and the angle 
ip is given by 

(p = argr(^^^) + ^ (^1 + log^^ . (44) 

The results in section VII of the paper correspond to setting A = in these formulas. 

Let us dehne the density of particles with dehnite momentum k, Uk = |/dfc(t)p. Using 
Eqs. (|38D and (^3]) and the normalization condition |a;(.p — = 1, we hnd the change 

in this density during the production event at tj: 


ni 


= e 


+ (1 + 2e--^)ni 

2g-7rA/2 + n{\fnl sin•d, 


where'd = (p + 291 + ~ cirga;(,. When the occupation numbers are large, 

be approximated by 


4 - 

where the “growth index” is given by 

1 


ni = exp(2 Vfc) ni 


hi 


= ^In (^1 + 2e“^^ - 2e“’^^/2-\/lTe^sinh) 


(45) 
]) can 

(46) 

(47) 


This growth index characterizes the “efficiency” of each particle production eventIt 
was argued in that in an expanding universe the phases '9 corresponding to different 
production moments tj are practically uncorrelated, and therefore it is a reasonable 
approximation to treat the phase as a random variable. The effective growth index pk 
can then be obtained by averaging over the phase h, leading to the expression 


nk{t) ^ -exp(2m/ifct). 


(48) 


It was demonstrated in that the results of this approach agree reasonably well with 
more exact numerical results. 


It is clear from (^71) that successful particle production requires vrA ^ 1. For generic 
initial conditions on the inflaton, 9 3> \Jm/M, this condition practically coincides with 
the condition (^ij) obtained from the more qualitative analysis of section [4.1| . For special 
initial conditions, 9 ^ ^Jm/M, the inflaton passes through the instability region and 
the quantity A becomes negative at sufficiently low momenta. Since the growth index 
increases with decreasing A, this leads to enhanced particle production rate. (In the 
non-supersymmetric case, the maximal production rate is achieved when the bare mass 
of the matter scalar is set to zero, corresponding to A = > 0. This is the case 

considered in [^.) For example, if sin6' = 0, the minimal value of A is about —mM/k‘1 ~ 

®Note that for some values of 9 the growth index could become negative, corresponding to a decrease 
of the particle number at time tj. 
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—M/\Jg^SQ — M^. If M <C gSo, the enhancement is insignihcant. For gSo ~ M, 
however, this effect can become very important. 

The rates of X 2 particle production can be derived in the same way. The only 
difference is the sign change in the last term of (0). As a result of this change, A 
is positive-definite in this case. At any rate, as long as M -C gSo and m -C M, the 
produced number densities of Xi and X 2 particles are nearly identical. 

Finally, let us note that the above analysis could also be performed in the non- 
super symmetric model of Section ||; one just has to use the corresponding expression 
for the effective frequency, Eq. (|]), instead of (|33D. The results are identical to the 
ones presented here with the replacement sin^ 9 — mM cos 9 —> Mq. (Of course this 
replacement does not make sense for the special initial conditions, 9 ^ ^Jm/M, in the 

supersymmetric model.) In particular, for generic initial conditions, 9 S> yjm/M, the 
production rates of Xi 2 bosons are the same as for a particle of the bare mass M sin 9 
in the non-supersymmetric model. 


A.2 6 7^0. 


Now, consider scalar preheating in the case e 7 ^ 0. As we discussed in Section ^ 3 , non- 
adiabatic particle production can only occur when the spiral inflaton trajectory passes 
through the instability region (plSf) . Let us consider a circular trajectory satisfying this 
property: Si = |5*|sinmf, S 2 = IF^Icosmt. The effective frequency for the held Xi, 
dehned in Eqs. (|3^) and is given by 


= 


a^{t) 


+ 2M^ (1 — sin mt) — mM, 


(49) 


where we have used S'* = —M/g. The inhaton passes through the center of the instability 
region at times tj = m“^( 7 r /2 + jn). Since the instability region is small, each of the 
particle production periods is short compared to the inhaton oscillation period m~^. 
This allows us to Taylor expand the sine function in for t K,tj\ 

k2 


= 


a\t) 


mM + M^m^ {t — tjY. 


(50) 


Dehning r = \/Mm{t — tj), the mode equation (|3^) at t = tj becomes identical to (^OD, 
with 

+ o?{tj)mM' 

Thus, the analysis of particle production in the e 7 ^ 0 case ehectively reduces to the 
analysis for e = 0 performed in the previous subsection. The change in the occupation 
numbers Hk during a single production event at time tj is given by Eq. (|^); assuming 
large Uk we obtain 


(1 -I- 2 e — 2 e sm'd)n{. 


(52) 
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where "i? is a random phase and A = A;[. This formula and Eq. ( |^ make it clear 
that particle production by an inflaton on a spiral trajectory is extremely efficient: the 
occupation numbers of states with low momenta typically grow by a factor of ~ 25 
in a single production event. A rough estimate of the energy density of the produced 
Xi particles yields 


pit) ~ 


m?M‘^ 




(53) 


where N is the number of production events before t. This estimate suggests that this 
process is not, in fact, accurately described by the stochastic resonance picture: after 
just two or three production events, the energy density of the X particles becomes 
comparable to that of the inflaton, and their backreaction becomes crucially important. 
Such explosive particle production in the instability region is similar to the tachyonic 
preheating which occurs in hybrid inflation models. 


B Ellipticity of the Inflaton Trajectory. 

When e = 0 in the superpotential ('|lT|), the complex phase of the inflaton held S is 
conserved, and its trajectory in the complex S plane is a straight line. When e is turned 
on, the phase is no longer conserved, and the trajectory is more complicated. As we have 
argued in section the evolution of the inhaton can be divided into two stages: during 
the hrst stage the phase-violating terms are important, while during the second stage 
they can be neglected. The inhaton trajectory during the second stage is an elliptical 
spiral; its ellipticity is determined by the amount of angular momentum accumulated in 
the hrst stage. In this appendix we would like to argue that even for small violations of 
phase invariance (e ~ 10 “®) the accumulated angular momentum is sufficient to obtain 
order-one ellipticity. 

The inhaton potential has the form 

V =\mS + e^ 2|2 _ ^ 2|^|2 ^ ^rn\S\‘^{S + S*) + e^l^- (54) 

The Noether current for the phase invariance is 

ns = J d^xi{S*S - S*S). (55) 

Both the hrst and the third terms in the potential (|5^) preserve the phase invariance, 
while the second term explicitly breaks it. In the presence of the expansion of the 
universe, the Noether current satishes the equation 

ris + 3Hns = t = temlSI^(S - S*). (56) 

In the case e = 0, the motion of the inhaton is described by 

S(t) = —So sin mt = — |S'o|e*^ sinmf, (57) 


22 





where 6 is determined by initial conditions. (Here we assumed that the motion of S 
has virialized, which is true when the amplitude is less than the reduced Planck scale.) 
For sufficiently small e, the violation of the phase conservation is small, and the inflaton 
trajectory is well approximated by (|57D . In this case, (^6]) can be integrated: 


poo 

nsa^{oo) = / dtiem\S\‘^{S — S*)a^{t) 

Jto 

roo 

= 2em\So\^altosm9 / dt-sm^mt, 

Jtn t 


(58) 


where a is the scale factor. The last integral is order unity. (If the integration is taken 
from to = 0, it is 7r/4.) 

The quantity determining the ellipticity of the inflaton orbit is 


ns 

mIS'P 


~ 2em|S'o|^aotosin6' 


m|S'pa^ 


2e|S'o|to sin^. 


(59) 


Because to is dehned to be the moment where S starts oscillating around the origin 
following the harmonic oscillator potential, to — {m\So\) (M* = Mp/\/^), we hnd 


ns 


~ 2e- sm 9. 

m 


(60) 


The above analysis is valid only if the ellipticity is small. For generic initial conditions 
on the inflaton and m ~ 10^^ GeV, this implies e -C 10“®. The analysis breaks down for 
e ~ 10“®, indicating that the ellipticity of the orbit in the case studied in Section 4^ is 
of order one unless the initial conditions are tuned so that 0 ~ 0. 
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